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In this paper we consider a class of semihamiltonian systems characterized 
by the existence of a special conservation law. The density and the current of 
this conservation law satisfy a second order system of PDEs which has a natural 
interpretation in the theory of flat bidifferential ideals. The class of systems we 
consider contains important well-known examples of semihamiltonian systems. 
Other examples, like genus 1 Whitham modulation equations for KdV, are related 
to this class by a reciprocal trasformation. 



1 Introduction 



Bidifferential ideals play an important role in the theory of finite-dimensional 
integrable systems, in particular in the bihamiltonian theory of separation of 
variables [20115]. 

Some recent results [IJ [19] suggest that they have also some applications in 
the theory of infinite-dimensional integrable systems, in particular in the case of 
integrable quasilinear PDEs. 

In this paper, following [19], we want to deepen the study of these applications 
in the case of diagonal integrable systems of quasilinear PDEs, the so called 
semihamiltonian systems. 
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Definition 1 \26f A diagonal system of PDEs of hydrodynamic type 

v i t = v l {u)u i x i=l,...,n (1) 

is called semihamiltonian if the coefficients v l (u) satisfy the system of equations 



M^)=Ml&) **J*k*i, (2) 



d k v% )=d k ( 9jV 

3 \v i — v k J \v l 

where di = -^ 



du l 

The equations (|2j) are the integrability conditions for the system 

djW 1 djV 1 

yjl — yjj yl — yj 



(3) 



which provides the characteristic velocities w l of the symmetries of ([I]), and for 
the system 

(V - v j )d t djH = d iV j djH - djv'diH, (4) 

which provides the densities H of conservation laws of (pQ). 

The knowledge of the symmetries of the system (UJ allows one to find its 
general solution. Indeed, according to a general scheme of integration of semi- 
hamiltonian systems proposed by Tsarev, the generalized hodograph method 
|26j . any solution of a semihamiltonian system is implicitly defined by a system 
of algebraic equations 

w l {u) = x + v l (u)t i = l,...,n (5) 

where the functions w l (u) are the solutions of the system ([3]) . 

A classical result in the theory of first order quasi-linear PDEs [18] states 
that, if system ([T]) possess a conservation law 

d t H + d x K = 

then the characteristic velocities v l can be written in the form 

d,K . 

i = l,...,n. (6) 



9iH 



This result has some interesting consequences in the case of a semihamiltonian 
systems. Due to integrability conditions ([2]) the space of solutions w 1 of the 
system ([3]) is parametrized by n arbitrary functions of one variable. 

Since the system @ is invariant with respect to the substitution v l — ► w l , for 
any solution (w , . . . , w n ) of the system ([3]) there exists a function K' such that 

In other words the characteristic velocities of the symmetries can be obtained 
applying the linear operator 

v h(') ■■= -^H 9 ^) i = l,-,n, (8) 
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to a suitable current K' . 

Note that, in terms of the density H and of the currents K and K', the system 
of algebraic equations © reads 

d(K' + xH - tK) = 0. 

Substituting ((7j) in (|3]) and taking into account (j4|) we obtain the equations for 
the currents: 

di d 3 K> = |^^^' - **_d^ 8K , m (g) 

In general, the problem of finding the solutions of the system ([9]) could be 
very difficult. The aim of the present paper is to study a special class of semi- 
hamiltonian systems characterized by the existence of a density of conservation 
law H such that the equations ([9]) for the associated currents reduce to the form 

(f - f j )didjK' = d ig %K' - djgidiK', (10) 

where /* = f l (u l ) and g % = g l {u l ). 

Surprinsingly, also the density H is a solution of the system (I10p . Therefore 
the solutions of (fTUj) play a double role: 

- fixed H, they are in one-to-one correspondence with the symmetries of a 
semihamiltonian system (see formula ©). In other words they define a semi- 
hamiltonian hierarchy 

- they label these hierarchies: different choices of H correspond to different 
hierarchies. 

The theory of flat bidifferential ideals arises naturally in this framework. First 
of all because any solution of the system (llOp defines a flat bidifferential ideal. 
Second because it provides a recursive procedure to compute the solutions of 

(P33J- 

The paper is organized as follows: in section 2 we recall some useful results 
about the theory of bidifferential ideals. In section 3 we apply these results to 
the theory of semihamiltonian systems. Section 4 is devoted to a discussion of 
the Hamiltonian formalism. In particular we find a class of metrics satisfying 
a system of Egoroff-Darboux type. Remarkably, in general, these metrics are 
not related to any Frobenius manifold, since their rotation coefficients are not 
symmetric. Finally, in section 5, we put reciprocal transformations into the game. 

Aknowledgements. I wish to thank B. Dubrovin, G. Falqui, F. Magri and 
M. Pedroni for useful discussions. I am also grateful to M. Pavlov for interesting 
comments. This work has been partially supported by the European Community 
through the FP6 Marie Curie RTN ENIGMA (Contract number MRTN-CT- 
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2 Bidifferential ideals 

A tensor field L : TM — ► TM, of type (1, 1) on a manifold M, of dimension n, 
is torsionless if the following identity 

[LX, LY] - L[LX, Y] - L[X, LY] + L 2 [X, Y] = 



is verified for any pair of vector fields X and Y on M. According to the theory of 
graded derivations of Frolicher-Nijenhuis [llj . a torsionless tensor field L of type 
(1,1) defines a differential operator d,L, of degree 1 and type d, on the Grassmann 
algebra of differential forms on M, verifying the fundamental conditions 

d-d L + d L -d = d 2 L = 0. 

On functions and 1-forms this derivation is defined by the following equations 

d L f(X) = df(LX) 

d L a{X, Y) = Lie LX (a{Y)) - Lie LY {a{X)) - a{[X, Y] L ), 



where 



[X, Y] L = [LX, Y] + [X, LY] - L[X, Y}. 



For instance, if L = diag(f 1 (u 1 ), . . . , f n {u n )), the action of dz on functions is 
given by the following formula: 

We can now define the concept of bidifferential ideal of forms. 

Definition 2 A bidifferential ideal 3 is an ideal of differential forms on M which 
is closed with respect to the action of both d and d^: 

d(3) C J d L (J) C 3 

For instance, if the ideal 3 is generated by a single 1-form a, the condition of 
closure with respect to the action of d and di reads 

da = A A a, d^a. = fj, A a, 

where A and \x are suitable 1-forms. 

In this paper we need a special subclass of bidifferential ideals, called flat 
bidifferential ideals. 

Definition 3 A flat bidifferential ideal 3, of rank 1, on a manifold M endowed 
with a torsionless tensor field L : TM — > TM, is the ideal of forms generated by 
the differential dh of a function h : M —>M obeying the condition 

ddLh = dhAda (11) 

with respect to a function a which satisfies the cohomological condition 

dd L a = 0. (12) 

Remark 1 In the language of Dimakis and Muller-Hoissen, the pair (d, di + 
da A) defines a "gauged bidifferential calculus". Some applications of this calculus 
to the theory of integrable systems are discussed in fl^. 



From now on, if not stated otherwise, we assume that the eigenvalues of L 
are pairwise distinct. In this case the general solution of the equation (|12p is 
given by the sum of n arbitrary functions of one variable: 



n 



£*'(* 



and the cohomological equation (jlip reads 

(f ~ f)didjh = djhdig' - dihd^gi (13) 

where f l {u l ) and g" 1 = g l {u l ) are arbitrary functions of one variable. 

3 Flat bidifferential ideals and semihamilto- 
nian systems 

In this section we show that the linear operator ([8]) establishes a one-to-one 
correspondence between the space of solutions of the cohomological equation 
(lllj) (which is parametrized by n arbitrary functions of one variable) and the 
space of symmetries of a semihamiltonian system. Indeed, it is easy to prove the 
following proposition. 

Proposition 1 Let H(u) be a solution of the cohomological equation OTjJ . then: 

1) the systems 

ut = v i H (K 1 )ui, i = l,...,n, (14) 

u\ = v l H (K 2 )u l x , i = l,...,n (15) 

commute for any pair (K\,K<i) of solutions of lill\) . 

2) the system of quasilinear PDEs 



u\ = [vh(K)K 



diK 



l,...,n (16) 



diH 

is semihamltonian for any solution K of the equation (EJ). 

Proof. 

1) The commutativity condition for the systems (|14p and (|15p reads: 



v}j(Ki) - <(^i) v l H (K2) ~ v> H (K 2 ) 

By straightforward computation we get: 

djv\j(K) ^djHdidjHdiK - d i Hd l d j K _ d ia djH 
v)j{K)-i? H {K) ~ diH djKdiH - diKdjH ~ p - /* flfefl" 

which does not depend on K. 



(17) 



2) By definition of semihamiltonian system we have to check that the character- 
istic velocities v l H {K) satisfy the system ([2]). For % ^ k ^ j ^ i, we obtain the 
identity: 

d f Ojvhm \ f *, 8 jH 



v\j{K) - v> H {K) ) \p-fdiH 



dj,a 



{d l Hf 



djadiHdkH dkadiHdjH diddjHdkH 



(p-p)(p-p) (r - n(p - n (p-p)(f k -p) 



which is clearly symmetric w.r.t. the indices j and k. 



□ 



We have constructed a family of semihamiltonian systems depending on 
functional parameters: the eigenvalues P(u l ) of L, the functions a and H. If 
difi 7^ 0, (i = 1, . . . , n), without loss of generality, we can assume p(u l ) = u % ; 
the two cases being simply related by the change of coordinates u l — ► P(u l ). 

Clearly, in order to make effective the construction one has to solve the coho- 
mological equation (jlip . Even if its general solution, depending on n arbitrary 
functions of one variable, is known explicitly only in some special cases (see sec- 
tion 6 of [22] and references therein), the double differential complex defined by 
the pair (d, d£) allows one to construct iteratively a countable set of solutions. 

Lemma 1 Let K$ be a solution of ([77]) . Then, the functions K[ defined recur- 
sively by 

dK l+1 = d L K { - K t da, (18) 

satisfy the equation ill]) . 

The proof is based on standard arguments in the theory of bidifferential ideals 
(20|, fTj . We report it for convenience of the reader. 

Let us start with the first step of the recursive procedure 

dK x = d L K - K da, (19) 



First of all, let us verify that the 1-form appearing in the right hand side of (I19j) 
is closed. Indeed, since Kq is a solution of (fTT]). applying to the right hand side 
of (I19H the differential d we obtain 



d (diKs - K da) = dK Q Ada- dK Ada = 0. 

So the function K\ is (locally) well defined. Moreover 

ddiK\ = dLKo Ada = dK\ A da. 

We prove now the theorem by induction. Suppose that 

dKi = d L Ki_ x - Ki_ x da 
dd^Ki = dK h A da. 
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Then the 1-form in the right hand side of (I18J1 is closed: 

d (d L Ki - Kida) = dK x A da - dK t A da = 
and satisfies the equation (fTTI) : 

ddiKi + \ = diKi Ada = dKj+i A da. 

□ 

Let us illustrate how to apply the previous procedure in the case H = a, 
Kq = —a. Using the recursive relations (fTHj) we get 

i d i R i i 

dta 

< = -^< = [f i -a}ui = [f i + K ]ui 

o^a 



o%a 
Following [TS] we can write the above hierarchy in the coordinate-free form: 



"to = 


= Ux 


Uti ~- 


= [L + K Efu x 


Ut2 -~ 


= [L 2 + K L + K x E]u x 



u tn = [L n + KqL"- 1 + K^' 2 + • • • + K n „ x E]u x 

where u is the column vector (u , . . . , u n ) t , E is the identity matrix and L is a 
torsionless tensor field of type (1,1). The above vector fields commute also in 
the non diagonalizable case |19j . 

Example 1 H = a, L = diag(u l , . . . ,u n ), a = cTr(L) 

Kq = —a = —c y^ u 3 

j 

\ 3 

* 2 = -f B*') 3 + 7 B*') 2 £«''-£(£•■ 

i i j 

and so on. 



Example 2 (Non diagonalizable case). Let 



L 



u 3 


a 2 
2 








U 3 


u 2 
2 








u 3 



The function a = u l {u 2 ) 2 satisfies the cohomological equation dd^a = 0. There- 
fore the first non trivial flow of the hierarchy starting from Kq = —a is 



u t 



II) 



a T 










2 







u 



U 2 



«' 



TTie other non trivial flows can be obtained solving the recursive relations [18] 
for the functions K\, K 2 , ■ ■ ■ : 



K 2 



-u\u 2 ?u 3 -\(u 2 f + \^?{u*f 

-u l {u 2 ) 2 (u 3 ) 2 - \{u 2 ) A u 3 + (u 1 ) V) V + 



and so on. 

Remark 2 Semihamiltonian systems of the form 



(f-a)u 



(20) 



have been obtained in [22] as finite component reduction of an infinite hydrody- 
namic chain. The connection between bidifferential ideals and such systems has 
been investigated in }19j . The starting point of that paper was the observation 
that the conditions OJ) and (dP for systems IWi) coincide with the cohomological 
equations 

dd L H = daAdH 
ddia = 0. 



4 Some remarks about the Hamiltonian struc- 
ture 

The Hamiltonian formalism for systems of hydrodynamic type was introduced by 
Dubrovin and Novikov in [21 0] . They considered first order differential operators 
of the form 

pv= 9 v(u)d x - g is r j sk (u) u k x (21) 

and the associated Poisson brackets 



1 J J 5u l 5ui 
where F = J g{u)dx and G = J g{u)dx are local functionals. 



(22) 



Theorem 1 J3}/ If det g l i ^ 0, then the formula |Hj) defines a Poisson bracket if 
and only if the tensor g lJ defines a flat pseudo-riemannian metric and the coef- 
ficients T J sk are the Christoffel symbols of the associated Levi- Civita connection. 

Non local extensions of the bracket (J22H . related to metrics of constant curvature, 
were considered by Ferapontov and Mokhov in [8]. Further generalizations were 
considered by Ferapontov in [7]. 

Let us focus our attention on semihamiltonan systems ([1]) 

u\ = v l (u)u x , i = 1, . . . n. 

In [7] Ferapontov conjectured that any semihamiltonian system is always Hamil- 
tonian with respect to suitable, possibly non local, Hamiltonian operators. More- 
over he proposed the following construction to define such Hamiltonian operators: 
1. Find the general solution of the system 

c\ i 

d 3 \n^g-=-^-. (23) 

v J — v % 



To this purpose is sufficient to find one solution ga of (j23f) . Indeed, the general 
solution is — rrr\, where (p l are arbitrary functions of one argument. The flat 
solutions of (|23p provide the local Hamitonian structures of the system ([1]) . 

2. Write the non vanishing components of the curvature tensor in terms of 
solutions w l a of the linear system ([3]) : 

4 = E e °« e « = ±1 - ( 24 ) 

a, 

(Ferapontov conjectured that it is always possible to find the expansion (|24|) ). 

Then the system ([1]) is automatically Hamiltonian with respect to the Hamil- 
tonian operator 

P*i = g^d x - g^uyl + £ e Q «d-VX. (25) 

a 

Example 3 |7| \27[ [21\j Let us consider the system of chromatography equations 
in Riemann invariants 

4= pjju*]" 1 *** i = l,...n. (26) 



The general solution of \23\) , in this case, is 



UkM u 



k _ ui) 2 



Qii = ~, — ^ I = i-, . . .71. 

y ^{u l ) 

where ip' l (u l ) are n arbitrary functions of one variables. For n > 3 all these 
metrics are not flat }21^ . They generate nonlocal Hamiltonian operators of the 
form [?]/ 



n 



pv = g «5Vd x - g is Ti k (u)u k x - Y, Kv^d-'wiui, (27) 



a=l 



where 



W: 



0; 



wt = a 



><£" 



Note that system &26\) can be written in the form 

diK 



u t 



~ dry « 



i = 1, . . . n, 



(28) 



(29) 



where 



K 



is a solution of the cohomological equation £77]) with L = diag{u l , . . . u n ) , a 

-Tr{L). 



Example 4 J2j/ Let us consider the semihamiltonian system 



Ut 



^y+2«* 



.4 = 1 



"X, 



. n, 



(30) 



T/ie general solution of \23\) is 

Uk^i(u k 



9i. 



^{u 1 ) 



1, 



n, 



where (f l (u l ) are n arbitrary functions of one variables. The choice tp l {u l ) = {u l ) a 
(a = 0,...,n) provides n + 1 flat metrics. For generic ip l {u l ) the metric ga is 
not flat and generates nonlocal Hamiltonian operator with infinite nonlocal tail. 
Note that system [W\) can be written in the form l2~9i) where 

2 



K 



iD»') 2 + i(E" 3 ) • 



is a solution of the cohomological equation ([77]) with L = diag{v} , . . . u n ) , a 
-\Tr{L). 



Let us consider semihamiltonian systems of the form (|16p . Taking into ac- 
count the equation (fTT]) . the system (f23l) reduces to: 



^-W-frr 



(31) 



Prom now on, in this section, we assume a = cTr(L) = c^™ =1 / J . In this case 
the general solution of (|3"Tj) is 



9i 



{d t Hf 



vwmvitf* - f l )\ 



lX\2c- 



l,...n, 



(32) 



where (p l {u l ) are n arbitrary functions of one variable. 

The rotation coefficients of the metrics (|32p depend on the constant c, the 
eigenvalues (J 1 , . . . , f n ) and on the choice of the arbitrary functions (y>i, . . . , <p n ) 
but not on the function H. More precisely we have the following proposition. 
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Proposition 2 Let H be a solution of the system: 

(f - f j )didjH = cdifdjH - cdjfdiH 



(33) 



Then, if a = cTr(L), the rotation coefficients flijiu) = — r jj - of the metrics 



( fffffj) . does not depend on H. More precisely they are given by the following 
expression 

Um(f-f l ) 



A- 



Ui^(P-f) 



'Vi 



P-f'Mfj 



Proof. 



(34) 



/%(«) 




rwp-/ z 



didjH djR djfj 



d { H diHfi-fi 



Taking into account (|33j) . we obtain formula (|34p . 



□ 



The problem of finding the expansion (|24p for the metrics (|32p is, in general, 
very difficult. 

For H = a this problem has been solved only for c = ±1, 1/2 (see section 9 
of [22J). 

The case H ^ a, which to our best knowledge, has not previously considered 
in the literature, can be reduced to the case H = a = cTr(L). Indeed, we have 
the following proposition. 

Proposition 3 Let H be a solution of the system 

(f ~ f j )didjH = cdifdjH - cdjfdiH 

and 



lit 



djK a • 



(35) 



the corresponding semihamiltonian hierarchy constructed with the solutions K a 
of the system (JHJ). Suppose that the hierarchy 



*t a — ^a^x 






is Hamiltonian w.r.t. the Hamiltonian operator 



(36) 



(37) 



and the Hamiltonian densities h a . Then the hierarchy i35\) is Hamiltonian w.r.t. 
the Hamiltonian operator 



pa = f 5 *3 dx _ §«fj fc («)«* - Y,« d - 



r. W a U xi 



(38) 
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where g % 



C d,H 



g n , the coefficients Y\ k are the Christoffel symbols of the 



associated Levi-Civita connection and w l a = C qh w l a . Moreover the Hamiltonian 
densities h Q of the systems 15*5)) can be obtained from the Hamiltonian densities 
h a solving the compatible system 

dih a = -?—dih a . (39) 

dia 

Proof: the non vanishing components of the curvature tensor 

can be written in the form 

E$ = - S ij (41) 

where the quantities S^ do not depend on H . Indeed in terms of the rotation 
coefficients (that do not depend on H), formula (p5)j) reads: 

4 = -4^4= ( di^ + d^ + x; fa?* 



Using this fact it is easy to obtain the expansion (|24l) for the non vanishing 
components of the curvature tensor of the metric g' u = c 2 ( -^-p ) gw- 

* diHdjH ^ 

Observe that the coefficients w l a = c-^jw l a = — J £ are characteristic veloc- 
ities of symmetries of (|35p . Therefore the bivector (|38p satisfies all Ferapontov 
conditions. Indeed: 

- the diagonal metric gu is a solution of the system (|23|) . 

- the coefficients T\ are, by definition, the Christoffel symbols of the associ- 
ated Levi-Civita connection. 

- the non local tail of (|38|) is constructed with the characteristic velocities w l Q 
appearing in the expansion of the non vanishing components of the curvature 
tensor. 

We have to show now that the function h a are hamiltonian densities. First of 
all we observe that they are well-defined. Indeed the compatibility of the system 
reads 

which is nothing but the system ^ for the densities of conservation law of the 
semihamiltonian hierarchy (|36|) . Moreover it is easy to check that if the functions 
h a are solutions of the system (|42p . then the functions h a are solutions of the 
system (j4|) for the densities of conservation laws of the semihamiltonian hierarchy 

dSSD- 

12 



□ 



We conclude this section mentioning an important property of the metrics 
32]) in the case L = diag(u , . . . , u n ) and ipi = 1 (i = 1, . . . ,n). 



Proposition 4 If f l {u l ) = u l and tfi = 1 (i = 1, . . . ,n) ; i/te rotation coefficients 
P^P satisfy the system 



dkPij = PikPkj ij= 3 + k 

^2 dk &j = ° l ^j 



Y^ u k d k Pij 



-A, 



i.i 



i + 3 



(43) 
(44) 

(45) 



Proof: the equations (|43p are automatically satisfied because they are equivalent 
to the conditions ([2]). 

Moreover, by straightforward computation we obtain 

2 , ^fcAj = fy'Aj + SiAj + 2^ PikPkj = 
k k^i,j 






c-l 



u-? — u l 



iw«* - «') iw«' - ul ) e^- ru- fe(« j - ^' 



u/ 



k^jiu^j,, 



Ui^-u l 



U ¥j (^-w 



l\2 



-+ 



e** iw«' - «') , n^-^rw^-^) 1 



Ui^j(u j -u 



IY2 



u> — u 1 



£ 



i -j 



1 1 

+ 



U k — U l Ul — U k 



U ¥j (u j -u l ) 

Y[l^{u l -u l ) 
Wl^-n l ) 

_ >p i i y- i y- / i 

Z — J iij — 11.™ £- — ^ ll z — ll_k Z — J \ 7//C _ 



v? — u 1 



- X 



u? — u k * — ' W — W 



kj^ij 



+ 



u k — u % u> — u h 
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and 






i.l 




D 

Remark 3 In general, the rotation coefficients \34ty are not symmetric. In the 
case of symmetric rotation coefficients, the equations [Jty , (44\) an d (45\ ) arise 
naturally in the framework of Frobenius manifolds f^j. 

As well-known the theory of Frobenius manifolds is related to the theory of 
isomonodromic deformations. Indeed, the equations $3\ ), (44\ ) an d fiJSty are 
equivalent to the system: 



d k V(u) = [V(u),[E k ,T]}, 



k = l, 



,n. 



(46) 



where 



(Ek)ij — 5ik$kj 

U := diag{u x , . . . ,u n ) 

T(u) := (/%) 

V(u) := («< - ^')/%, 

that governs the monodromy preserving deformations of the operator 



d_ 

dz 



U + 



The proof of this fact (see /§[/,) does not rely on the hypothesis of symmetry of the 
rotation coefficients. 
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5 A remark on Whitham equations 

A well-known example of semihamiltonian system is the system of quasilinear 
PDEs that describes the slow modulations of g-gap solutions of the KdV hierar- 
chy: the Whitham equations [28J. In this case the characteristic velocities can 
be written in terms of hyperelliptic integrals of genus g. In g = 1 case these 
equations read: 



u 



1 + u 2 + u 3 2{u 2 -u l )K{s) 



+ 



u t 



U 1 + V? + u 3 



+ 



3(K(s)-E(s 
2(u 2 - u 1 )K(s 



a , 



3 ' 3(E(s)-(l-s 2 )K(s) 

u 1 +u 2 +u 3 2(u 3 - n 1 )(l - s 2 )E(s] 



n , 



a ; 



3E(s) 



where s = " 3 ~"i , K(s) and E(s) are complete elliptic integrals of the first and 
second kind. 

For Whitham equations the hodograph method is effective [15l [TBI E31 I2H 
25l [T3l [T^t El [12]. Indeed it is possible to construct explicitly the symmetries 
appearing in the equations ([5]): 



Theorem 2 J721 03J [^ E3 \T$ There exist functions qi(u),q2(u), q${u) such 
that the characteristic velocities w l of the symmetries of the Whitham equations 
have the form 

w i := [1 + qidi] K, i = l,..,3, (47) 

where the function K is a solution of the following system of Euler-Poisson- 
Darboux type: 



2(n* - u j )didjK = d { K - djK 
which can be explicitly solved. 



i ¥=3, h3 = 1,2,3. 



(48) 



The functions qi(u) can be written in terms of the complete elliptic integral K(s) 
and E(s). Moreover from the conservation of waves it follows that [17 |, I13 |, IT^ ]. 125] : 



H 

~i\H 



1,2,3, 



where 



d( 



H = & g (49) 

is the wavelength (the integration is taken over the cycle around the gap u 2 < 
£, < u 3 )- Therefore, the Whitham equations can be written in the form: 



v l H (K) 



K 



H 



diK 



1,2,3. 



Note that the wavelength H satisfies the Euler-Darboux-Poisson system (J48J), 
which is a particular case of the cohomological equation (jlip corresponding to 
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the choice L = diagiu 1 , u 2 , u 3 ), a = — ^Tr(L). This remark suggests to consider 
systems of the form 



4 = »Mif)4 



K ~m diK 



U X1 



l,...,n. 



(50) 



where the function H and K are solutions of the cohomological equations (1331) . 
It is easy to prove that such systems are semihamiltonian and that, fixed H, the 
systems 



4 = v % H {K x )u % x , 



i = 1, ...,n 
i = l,...,re, 



commute for any pair (Ki, K%) of ([33]) . 

This fact can be proved by straightforward calculation or simply observing 
that systems of the form (|50h can be obtained from the systems studied in this 
paper by means of a reciprocal transformation. 



Proposition 5 Systems 



% 



K ~m a ' K 



i = 1, ...,n. 



are related to the systems 



■ii t 



djK 
diH 



u, 



by the reciprocal transformation 



dx = Hdx — Kdt 
dt = dt 

The proof is a trivial computation. 

6 Conclusions 

In this paper we studied some applications of the theory of flat bidifferential 
ideals to semihamiltonian systems of quasilinear PDEs. 

The starting point of the present paper was the observation that for any 
semihamiltonian system there exists a linear differential operator that, acting on 
a suitable domain, provides all the symmetries of the system. 

We showed that for a special class of semihamiltonian systems this operator 
and its domain are completely characterized by the solutions of a cohomological 
equation. 

Moreover the theory of flat bidifferential ideals naturally provides a recursive 
procedure to compute the solutions of this equation. 
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